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The dual superconductivity is believed to be a promising mechanism for quark confinement. Indeed, what this 
picture is true has been confirmed in the maximal Abelian (MA) gauge. However, it is not yet confirmed in 
any other gauge and the MA gauge explicitly breaks color symmetry. To remedy this defect, we propose to use 
our compact formulation of a non-linear change of variables on a lattice. This formulation has succeeded to 
■ extract the magnetic monopole with integer-valued magnetic charge in the gauge-invariant way. In this talk, we 

present measurements of various correlation functions for the operators constructed from the CFN variables 
in SU(2) Yang-Mills theory. Some of our results reproduce previous results obtained in MA gauge, e.g., 
DeGrant-Toussaint monopole, infrared Abelian dominance and off-diagonal gluon mass generation. These 
studies preserve color symmetry, in sharp contrast to the conventional MA gauge. We argue the gauge fixing 
independence of these results and the implications to quark confinement. 
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1. Introduction 

Quark confinement is still an unsolved and challenging problem in theoretical particle physics. The dual 
superconductivity]^] is believed to be a promising mechanism for the vacuum of the non- Abelian gauge theory [jl]]. 
Indeed, the relevant data supporting the validity of this picture have been accumulated by numerical simulations 
especially since 1990 and some of the theoretical predictions ||] have been confirmed by these investigations: 
infrared Abelian dominance [§], magnetic monopole dominance [^] and non- vanishing off-diagonal gluon mass 
[0] in the Maximal Abelian gauge [g], which are the most characteristic features for the dual superconductivity. 
However, they are not yet confirmed in any other gauge and the MA gauge explicitly breaks color symmetry. To 
establish this picture in gauge invariant way, we need to answer how to define and extract the "Abelian part" 
from the original non-Abelian gauge field g/^ which is responsible for the area decay law of the Wilson loop 
average. The conventional Abelian projection [||] is too naive to realize this requirement. At the same time, 
we must answer why the remaining part in the non-Abelian gauge field srf^ decouple in the low-energy (or 
long-distance) regime. 

We propose to use a non-linear change of variables (NLCV) which was called the Cho-Faddeev-Niemi 
(CFN) decomposition 11 , 12, 13, 14] to remedy the defect of ordinary approaches. [ |15[ [TJ, [18|] We introduce 



a compact representation of NLCV on a lattice. The naive decomposition presented at the last conference was 
improved to extract the magnetic monopole with integer- valued magnetic charge in the gauge-invariant way. Jl9| ] 
Some of our results reproduce previous results obtained in MA gauge, e.g., DeGrant-Toussaint monopole, in- 
frared Abelian dominance. 

2. Lattice CFN variables or NLCV on a lattice 

We propose a formulation of NLCV on a lattice It is a minimum requirement that such a lattice formulation 



must reproduce the continuum counterparts in the naive continuum limit. In the continuum formulation [12, 10], 



a color vector field n(x) = (n^(x)) (A = 1,2,3) is introduced as a three-dimensional unit vector field. In what 
follows, we use the boldface to express the Lie-algebra ra(2)-valued field, e.g., n(x) := ha(x)Ta, Ta = \®a with 
Pauli matrices Oa (A = 1,2,3). Then the ra(2)-valued gluon field (gauge potential) A^(x) is decomposed into 
two parts: 

A M (x)=V M (x)+X M (x), (2.1) 
in such a way that the color vector field n(x) is covariant constant in the background field V^(x): 

= %[V]n(x) := a M n(x) - ig[V M (x),n(x)], (2.2) 

and that the remaining field X^(x) is perpendicular to n(x): 

n(x) -X M (x) = 2tr(n(x)X jU (x)) = 0. (2.3) 

Here we have adopted the normalization \x(TaTb) = \8ab- Both n(x) and A„(x) are Hermitian fields. This is 



also the case for (x) and X^ (x) . By solving the defining equation ( 2.2 ), the (x) and the X^ (x) . are obtained 
in the form: 

V M (x) =vl(x)+V^(x) =c ll {x)n{x)-ig-\d ll n{x),n{x)l (2.4) 
X M (x) = -ig- 1 [n(x),%[A]n(x)] (2.5) 

where the second term V^(x) := — ig~ l [<9^n(x),n(x)] = g^ 1 (d ll n(x) x h{x))aTa is perpendicular to n(x), i.e., 
n(x) ■ (x) = 2tr(n(x)V^ (x)) = 0. Here it should be remarked that the parallel part v]i (x) = c^, (x)n(x), (x) = 
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tr(n(x)A jU (x)) proportional to n(jc) can not be determined uniquely only from the defining equation (2.2), and 
the perpendicular condition of ( j2j| ) determines vl (x) and remainder part (x) . 

On a lattice, on the other hand, we introduce the site variable n x = h^Oa in addition to the original link 
variable U xfl which is related to the gauge potential Ay^: 



U x .n =exp(-j££ 



(2.6) 



where = (x + ju/2,/x) stand for the midpoint of the link 1 . 

In what follows, we use the blackboard boldface to express the field determined by the link variable. Note 
that n v is Hermitian, nj = n^, and U Xj ^ is unitary, U xfl = U~^. The link variable U X) ^ and the site variable n x 
transform under the gauge transformation II [10] as 

-> ei x u X)ll al +ll = , n x -> n x n x nl = ri x . (2.7) 

Suppose we have obtained a "link variable" V x ^ and X Xj ^ as a group element of G = SU (2) through 

7 XlM =exp(-iegVy Jt ), (2.8) 



'x,HJ> 

X X)fl =exp(-iegX Xili ). 



(2.9) 



Vy^are related to the ,™(2)-valued background field Vy » where V^n is to be identified with the continuum 
variable (2.4) and hence V XjiU must be unitary V Xtfl = V x }- A lattice version of defining equation (22) and (2.3) 
are given by 



Dj? [V]n, := e- 1 [V x ^n x+ll - n x V x>fl ] = 0, 
tr(n x X x ^)=0. 



(2.10) 
(2.11) 



The defining equation ( |2.10[ ) needs a lattice covariant derivative for an adjoint field. We adopt the midpoint eval- 
uation of the difference d^n x = £ _1 [n A+iU — n r ] = d u n r /+ G{z 2 \ therefore the continuum covariant derivative 
for the adjoint field up to &{e 2 ) at midpoint: 



.2 



e' 1 [V x ^n x+l j,-n x V x , IJ .] = d^n x > -/g[V. v ^,n v /] - — {gV. v ^,^n y } + G{z 2 ). 
The derivative ( 2. lOj ) obeys the correct transformation property, i.e., the adjoint rotation on a lattice: 

o( £) [V]n^Q I (4 e) [VK)fl^> 
provided that the link variable V x ^ transforms in the same way as the original link variable U x ^ : 



V' 



(2.12) 



This is required from the transformation property of the continuum variable V^(x), 3 see [10]. Therefore, we 
obtain the desired condition between n v and V Xj n . 



n x V x ^ — V X) (iHx+n- 



(2.13) 



In general, the argument of the exponential in (2.6) is the line integral of a gauge potential along a link from itoi + /j. We 
adopt this convention to obtain the naive continume limit of ff(e 2 ). Note also that we define a color vector field n(x) := n A (x)T A in the 
continuum, while n A := i^Ga on the lattice for convenience. 



2 The term '-j {gV^n, <?^n t /} is of the order ff(e 2 ), since „ in contimume limit is obtained as eq(2.4) and d^ny • ny =0 + ff(e). 
3 This indicates that V x ,)i in ( |2.8| ) is considerd as the link variable whose argument of the exponential is the line integral of a gauge 
potential along a link from x to x + fl . 
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The defining equation ( 2.13 ) for the link variable V x ^ is form-invariant under the gauge transformation II, i.e., 

n x"x,fi V x,ii"x+H' 

A lattice version of the orthogonality equation ( |2.3| ) given by equation( [2.1 1[ ) or 

tr(n x exp{-/£gX^}) = tr(n,{l - iegX^}) + 0{s 2 ) = + ff(e 2 ). (2.14) 

This implies that the trace vanishes up to first order of e apart from the second order term. Note that X x>ll is 
defined on the lattice site and transforms in the same way as n v : 



X 



X' 

x,n 



(2.15) 



so that orthogonality condition < J2- 1 1[ ) is gauge invariant. 



Then, we proceed to solve the defining equation ( [2.1 3\ > for the link variable V x ^ and equation ( |2.11[ ) for 
the variable X x ^, and express it in terms of the site variable n v and the original link variable U x> ^, as is the 
case that the continuum variable Vm(x) and X^(je) are expressed in terms of n(x) and A jU (a). Remembering the 
&x,n ~ ^x,n , Xx,u can be defined using link variables V x ^ and U x u contacting to site x, and linear 



relation '. 



combination of Vx_n nUx-ii,ii an d U Xjfl V Xtfl , are candidate to satisfy the required transformation property ( 2.15 ); 



= exp(-iegX x ^) 



(2.16) 



2 2 



where the relation for matrices exp(eA) exp(efi) = exp(eA + eB + e 2 [A,B]/2 + ^(e 3 )) and its inverted version 
of exp(eC + £ 2 D) = exp(£C)exp(e 2 D) + &(e ) are used. The parameter X = <p is selected so that X Xyil is 
determined to coincide with continuum expression up to ^(e 3 ). 



As for V X) n , on the other hand, the equation ( J2. 13[ ) is a matrix equation and it is rather difficult to obtain the 
general solution. Therefore, we adopt an ansatz (up to quadratic in n): 

V x ,n = U x # + an x U x ^ + fiUx^nx+n + yn x U x ^n x+ll , (2. 17) 

which enjoys the correct transformation property, the adjoint rotation ( 2.12| ). It turns out that this ansatz satisfy 



the defining equation ( |2.13| ), if and only if the numerical coefficients a,f5 and y are chosen to be y = 1 and 



a = p. Then, substituting the ansatz ( 2.17 ) with a still undetermined parameter a into equation ( 2.16 ), we 
obtain a = O+0(s 2 ) (see @)- 



Thus we have determined V x ,n and X X n up to an overall normalization 

V x ,fi = y x ,\i [U,n] = U Xtll + n x U Xtll n x+ll , 
X x ,n =X x ^[U,n] = V x t^U x ^ + U x ^ 



x,ii 



The unitary link variable V x ^ [U,n] and X X)fl [U,n] can be obtained after the normalization: 



9 Xifl [U,n] := V^/yj ^tr[V x \^}, X^ [U,n] := X^/yj -tr^X^]. 
3. Numerical simulations and generation of configuration of NLCV 



(2.18) 



We generate configurations of link variables {U X n} using standard Wilson action. The numerical simulation 
are performed on 24 4 lattice at j8 = 2.3, 2.4, 2.5 by thermalizing 15000 sweeps, and on 36 4 lattice at j8 = 2.5, 
2.6, 2.7 by thermalizing 18000 sweeps. 200 configurations are obtained every 300 sweeps. 
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M-YM 



Yang-Mills 
theory 



_St/(2)„x[«/(2)/t/(l)] 



M-YM 



Type I SU(2) m 



V < j3quivalent ^> V +nMAG 



Type II SU(2)„ rf 




YM+LLG 

global St/ (2) 



/ IVI-TIVI 

< ^equivalent^ > VtLLG+nMAG 



, gauge orbit n x — B K 
(CFN variable) 




Figure 1: (left)The relationship between the master- YM theory and the original YM theory, (right) NLCV via gauge 
transfroamtion. 



The NLCV on a lattice is obtained according to the method of the previous paper JT^]. Figure [l] shows 
the extended gauge symmetry in the master- YM for NLCV (left panel) and NLCV of SU(2) link variables via 
gauge transformations (right panel). The configuration of the link variable U x ^ and the color vector field 
has an extended gauge symmetry SU (2) m x [SU (2) /U (l)] e . The equivalent theory to the original YM theory is 
obtained by the gauge fixing which we call the new Maximal Abelian gauge (nMAG). We define a functional 
written in terms of the gauge (link) variable U x ^ and the color (site) variable n v ; F„mag\U \n;£2, 0] = tr(l — 
n jr t2 ?7 v ^ n A:+ ^ a L^ lU ), where we have introduced the enlarged gauge transformation: n U x ^ := £l x V x ^ n| + „ 
for the link variable U x ^ and ®n x := © x n x Q ' l ®l for an initial site variable n x . The gauge group elements Q. x 
and & x are independent SU(2) matrices on a site x. After imposing the nMAG, the theory still has the local 
gauge symmetry SU(2)?~f,, since the "diagonal" gauge transformation (0 = 6 does not change the value of 
the functional F„mag\U ,n;Q,0]. Therefore, the configuration of n x can not be determined at this stage. In 
order to determine n v , we need to impose another gauge fixing or a choice of the gauge of link variable U x ^ 
for fixing SU {2)^. The desired color vector field is constructed from the interpolating gauge transformation 
matrix V by choosing the initial value n x = 03 and n x := X C30] = n^a A , n x = tr[cA0 v CT30|]/2 (A = 1,2,3) 
where {& x } are given by gauge transformations that satisfy U Xifl = & x U^n For example, we choose the 

conventional Lorentz-Landau gauge or Lattice Landau gauge (LLG) for this purpose. The LLG can be imposed 
by minimizing the function Fllg[U;£1] = Y,x,n ~~ a U X:li ) with respect to the gauge transformation Cl x for the 
given link configurations {U X n}. 



4. Infrared Abelian Dominance and Mass generation of the off-diagonal gluon 



Using new variables through NLCV, we are now ready to study characteristic features of the YM theory 
for any choice of gauge fixing such as infrared Abelian dominance, magnetic monopole dominance and the 
non-vanishing off-diagonal gluon mass. 4 Our proposed decomposition extract the "Abelian part" V x n in any 
gauge fixing preserving the color symmetry. The conventional MAG fixed theory is reproduced as a special case 
of our formulation base on NLCV. To study the infrared Abelian dominance and the non- vanishing off-diagonal 
gluon mass in LLG other than MAG, the correlation function of the decomposed variable V x ; ijU and X x u has 
been measured. Left panel of figure || shows propagators Daa(x — y) = (A Xjli A y ^\, Dyv(x — y) = (V A%iU Vy )jU ) 

and Dxx( X ~ y) ~ ^ Y ^ Tl, ° r, °" r, ° nntontlolo o.-o Hofino/I oc llnV uonoHoc A . — _d_ A _ A t 



The gauge potentials are defined as link variables Ay ijU = ^§ Ax,n — A x>fX 



X,[l 



representation, 



2g£ 



. On the other hand, we can define the X Xj ^ in two ways, one is extracted from compact 

yt " 



and the other is from definition of the decomposition (2.1), 



The magnetic monopole dominance has been found using integer valued and gauge invariant magnetic monopole defined by our 
NLCV. This fact has been reported in lattice2006 H . 
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1 ■ a 



D24[VV]b2.3 m-i 
D24[VV]b2.4 
D24[VV]b2.5 h«-i 
D36[VV]b2.5 heh 

D36[VV]b2.6 

D36[XX]b2.5 h«h 

D36[XX]b2.6 mm 
D36[XX]b2.7 hh 
D24[XX]b2.3 hth 
D24[XX]b2.4 i-fH 
D24[XX]b2.5 



r (sqrtfa)) r (sqrtfa)) 

Figure 2: (left) correlation functions (O(x)O(y)) in the logalithmic scale, (right) rescaled correratin function 

In (r 3 / 2 G^(r;M)) 



&*?,n ~ y» / ,ti- Plotting of two types of Dxx(* — y) overlap for several lattice spacings (several j8s) , and the 
extracttion of the variable is consistent (see left panel of figure §). On the other hand, D AA {x — y) and D vv (x—y) 
overlap, and Dxx{x — y) is dumped more quickly for infrared region than Dyy(x — y). This implies that the 
infrared Abelian dominance is found in the LLC 

Next we study the mass of the decomposed fields from the correlation functions. The inverse Fourier 
transformation of the massive gauge boson propagator should behave for large r = \x— y\ as follows, 

So the scaled propagator r'/ 2 G^^{r,M) is proportional to e~ Mr , that is, the mass of gauge potential M, is 
obtained as the dumping factor of r 3 / 2 G^(r;M). In other words, the gradient of the linear fitting in the r vs 
ln(r 3 / 2 G MA1 (r;M)) plot gives the mass M. Right panel of figure || shows the plots of the scaled propagator of 
n and Vy,u. The distance r is measured in the unit of square root of the string tension y/Osr (= 440 MeV), 
and vertical axis is scaled in the logarithm to measure the dumping factor by the linear fitting. To determine the 
physical scale, the relation between j3 and lattice spacing e is obtained from [[H]]. The dumping of propagator of 
Xjc^ gives the mass Mx ~ 1.18 GeV, and the "Abelian part" Y x ^ indicates My ^ 0.48 GeV. These are consistent 
with study in MAG. [{J 



5. Summary and discussion 

We have proposed a new formulation of the lattice Yang-Mills theory based on the NLCV which was once 
called the CFN decomposition. This resolves all drawbacks of the previous formulation of the decomposition 
on a lattice. This compact formulation enables us to guarantee the magnetic charge quantization in the gauge 
invariant way and to extract the "Abelian part" and the "off-diagonal part" preserving color symmetry in any 
choice of gauge of the original YM theory. These features are sharp contrast to the conventional MA gauge 
and these studies. We have measured the correlation function (propagator in real space) in LLG. The Infrared 
Abelian dominance and the gluon mass generation have been found. These results are consistent with study in 
MA gauge. 
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